NASA/TM— 2004-212970 



Generalization of Equivalent Crystal Theory 
to Include Angular Dependence 

John Ferrante 

Glenn Research Center, Cleveland, Ohio 
Fredy R. Zypman 

Yeshiva University, New York, New York 


March 2004 


The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:/ / www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 


NASA/TM— 2004-212970 



Generalization of Equivalent Crystal Theory 
to Include Angular Dependence 

John Ferrante 

Glenn Research Center, Cleveland, Ohio 
Fredy R. Zypman 

Yeshiva University, New York, New York 


National Aeronautics and 
Space Administration 


Glenn Research Center 


March 2004 


Acknowledgments 


We would like to thank the Tribology and Surface Science Branch, NASA Glenn Research Center, 
where part of the work was performed. We also thank Dr. Phillip Abel for fruitful comments. 


This work was sponsored by the Low Emissions Alternative 
Power Project of the Vehicle Systems Program at the 
NASA Glenn Research Center. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / /gltrs. grc.nasa.gov 


Generalization of Equivalent Crystal Theory 
to Include Angular Dependence 


John Ferrante 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Fredy R. Zypman 
Yeshiva University 
Department of Physics 
New York, New York 10033-3200 


Abstract 

In the original Equivalent Crystal Theory, each atomic site in the real crystal is assigned an 
equivalent lattice constant, in general different from the ground state one. This parameter 
corresponds to a local compression or expansion of the lattice. The basic method considers these 
volumetric transformations and, in addition, introduces the possibility that the reference lattice is 
anisotropically distorted. These distortions however, were introduced ad-hoc. In this work, we 
generalize the original Equivalent Crystal Theory by systematically introducing site-dependent 
directional distortions of the lattice, whose corresponding distortions account for the dependence 
of the energy on anisotropic local density variations. This is done in the spirit of the original 
framework, but including a gradient term in the density. This approach is introduced to correct a 
deficiency in the original Equivalent Crystal Theory and other semiempirical methods in 
quantitatively obtaining the correct ratios of the surface energies of low index planes of cubic 
metals-(lOO), (110), (111). We develop here the basic framework, and apply it to the calculation 
of Fe (110) and Fe (111) surface energy formation. The results, compared with first principles 
calculations, show an improvement over previous semiempirical approaches. 


Introduction 

Equivalent Crystal Theory (ECT) is a method to calculate energies of atomic conglomerates. It 
has been very successful in predicting formation energies of bulk defects. Bozzolo, Ferrante and 
Smith (BFS) have extended it to multicomponent metals, and have successfully used it to predict 
the ground state structures and surface properties of various alloys. 

The basic idea of the original Equivalent Crystal Theory (OECT) rests on the fact that the total 
ground state energy of a quantum system is a functional of the electron density [1], Therefore, 
OECT requires that, at the atomic site under consideration, the density be the same for both the 
real lattice and the equivalent, defect-free lattice. That provides a value for the equivalent lattice 
parameter a eq at that site. The energy contribution from that site is obtained by using the 
Universal Binding Energy Relation (UBER) curve. The UBER provides binding energy as a 
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function of average atomic separation (or in this case, lattice parameter). In this work, we 
generalize OECT without losing its original flavor, that is, by keeping the assumption that the 
electron density must remain the same for both the real and the equivalent crystal. However, here 
we extend it to include non-local effects by requiring that both, the real and equivalent crystal 
densities be the same not only at the center of the atomic site into consideration, but also in a 
neighborhood around that point. Concretely, we will impose that the density at an atomic site and 
its gradient be the same for the real and the equivalent crystal. These density variations in the 
neighborhood can be estimated from recent ab initio calculation for surface binding curves for 
Fe. The reference crystal now includes both isotropic distortions and separation of the solid in 
two pieces nonnal to a given surface plane. 

In this new framework, the reference crystal at an atomic site will be assigned two (instead of 
one) parameters: the equivalent lattice parameter and the relative separation, y, between two 
hypothetical surfaces. The direction of the separation, and the separation itself, do not have to be 
the same for all the real atoms in the crystal. We will call that, potentially variable, direction z. 
The two parameters are obtained from two requirements, which are: the site density and its 
gradient at each location are assumed the same for both the real and the equivalent crystal. 


Development of the Ideas 

In this modified version of ECT (MECT), the reference crystal is not just a perfect lattice with 
a -t- ao, but, in addition, it incorporates a pair of surface-tenninated bulks with a ^ ao, separated by 
a distance yd, where d is the interplanar distance normal to a particular plane. 






a 


ydH 



Figure 1 . — Reference crystal. 


NASA/TM- 


- 2004-212970 


2 



Figure 1 shows the reference crystal in a generic configuration. Clearly, as opposed to the 
original ECT (OECT), not all atoms in this reference crystal contribute the same amount to the 
total energy. We will consider the local energy contribution coming from a reference atom in the 
(quasi) surface (QS), shown in red in the figure. The quasi-surface becomes a bona fide surface 
in the limit y — > + «, An atom in QS (AQS) contribution to the total energy is now a function of 
both the lattice parameter a, and the relative distance y, between the two QS. This differs in form 
from OECT which expands the energy in a series containing a volume term, plus two-, three-, 
and four-body terms. In this formulation, the next term beyond the volume’s would be the 
gradient of the density, and the next beyond that would be the curvature. 

In OECT, a was found by requiring that the total electron density at the center of a given atom be 
the same for both the real defect crystal and the reference crystal. That condition will remain 
valid in MECT, but we will impose a second constraint. The motivation for the additional 
equation is to also look at the density close to the center of the given atom, not only at the center. 
We will impose that the real and reference densities must be equal in that small (but not zero) 
region. As in OECT, we take the densities from each atom to be given by 

Z | — * P — cx\ R — r\ • • • • 

\R-r e , with standard notation. Consider the atom for which we want to 

R 

calculate the energy contribution, to be located at the origin. We expand R-r taking into 
account that r « R: 

p(r ) = p(0) + r ■ V p(r)|_ =5 + • • • (1) 

Then, expanding the power and the exponential, 

p(F) « Y, RPe ~ aR +? X iaR ~ P) RP ~ le ~ aR I (2) 

R R 

From this expression, we will impose: 

Y J R p e~ aR = Y 4 R p e~‘* ( 3 ) 

R equivalent R real 

•n e = Y J ( aR -p) RP ^ e ^ aR i -K ( 4 ) 

Rreal 

where h e and h r are fixed directions in the respective spaces. 

Both right-hand sides are explicit functions of y and a and will be obtained by solving ( 3 ) and 

( 4 ). 


'£(aR-p)R p - l e‘*i 

R equivalent 
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Generalization of UBER 


The most critical issue to stress in MECT is the fact that, as opposed to what happens in OECT, 
not all atoms in the reference crystal are equivalent and therefore contribute a different amount to 
the total energy. 

Let us assign the energy contribution from each atom to the total energy 

E,=E Aes +(E B -E Aes )^-e- u -»\ (5) 

where all atoms in the same plane n will be weighted equally, but atoms in different layers will 
have a different contribution. Here E B is the contribution from bulk (which is given by the OECT 
UBER), Eaqs is the (still unknown) contribution from AQS, X ~ 1/M (M is the number of layers 
in which atoms are substantially different from bulk), finally n is the layer’s depth (n = 1 
corresponds to the top layer, our QS). Here we have assumed an exponential form for the 
contribution to the energy as a function of depth. Hung et al [2] have shown that the electron 
density in slabs of increasing widths rapidly approaches its bulk value . 

The total energy E tot can be calculated as: 


E 


Tot 


n, E E „ = N *- N c E °+ N ,' ( E w - E . ) • E 

n = 1 n = 1 


(6) 


Where N A is the number of atoms in-plane, N D is the number of layers (although the sum goes to 
infinity here, in practice Nd is very large but not infinite). 

Calling N = N a .Nd, the total number of atoms, we obtain: 


E Tot ~ N ' E B + 


N A ' ( E AQS E b) 


-X 


( 7 ) 


Hung et al [3] have done a first-principles density functional theory calculation where they 
obtained the energy as a function of separation of two surface-terminated bulks. Call that energy 
curve E(ydo), where do = cao is the equilibrium interplanar distance, ao is the equilibrium lattice 
parameter, and c is a geometric factor that depends on the orientation and the crystal structure. 
Ferrante and Smith [4] defined the adhesive energy as the energy two surfaces have referred to 
the energy they have when infinitely separated. Since reference [1] reports energy per unit area 
(A), we have, for one QS: 




( 8 ) 


We now make the identification E To t = 2E in equations (7) and (8), since they both correspond to 
total energies of the two QSs, calculated independently. Then, 
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p ( v \ = 2 N a E AQ s(Y) E/s Eaqs(+°°) E b 
ad(7 A [ \-e~ x l-e~ A \ ( 

When y = 1, the two QSs are separated by their equilibrium distance, and then X = 0 and 
E A qs = E b . Thus, in that case: 

2 N a E R -E An J+°°) 

E «m= / (io) 


Substituting (11) into (9): 


E. QS (+°°) = E B +A-(e- l -X)E ad (\) 


2 N , E, n Ay)-E„ 

E,Ar)=—r 

A \ — e 


Solving for E A qs(y), which is the energy of the reference atom AQS, 

E tQS (y) = E, + ^ f (\-e t )VEM)-E.m 


Call T the area occupied by an AQS, that is T = . If E B is the energy per atom, then the 


energy per 2D Wigner-Seitz area is: 


EjnAr) =^-+W- e~ x )[E ad (r) - E ad ( 1)] 


The tenn E B has been parameterized [5] by the Universal Binding Energy Relation (UBER): 


E D (x) = E 1 - 1 +a e 

B v 7 c 


where x = — , is the relative lattice parameter of the reference crystal, E c is the bulk cohesive 

a 0 

energy, and: 

a ' = 7 L (to)‘ (x “ 1) <16) 

B 

for BCC structures. Here is related with the bulk modulus and is known for metals [6,7], and a 
is the lattice parameter. 


E ad (y) in equation (12) has also been parameterized by a UBER in reference [1]: 
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(17) 


E,Ar)=-2n i+f (r-0 

It is important to note that the last expression is valid only for a = ao since that is the calculation 
performed in reference [1]. We need to extend it to any value of y and a. 



If they had done the same (adhesion) calculation for a crystal with a larger lattice parameter, the 
energy (for a given relative separation) would have been smaller. In particular, in the limiting 
case a — > + oo, the energy of the system tends to zero, for any separation. On the other hand, the 
first principles calculation A(y) (explicitly written in tenns of the relative separation) has a 
minimum for y= 1, that is for a separation equal to the lattice parameter (or the corresponding 
interplane distance). Finally, we write the excess (referred to the bulk ground state) energy of the 
reference atom AQS as: 


-fi-V**] +y s \ 1- 1 + -7^0-1) e 


This expression is built in such a way that the ground state of the system corresponds to y= 1 
and x = 1, that is an undistorted crystal with the equilibrium lattice parameter. Any 
configurational change, away from that structure, can only increase the energy. 


Setting the Reference QS (100) 

In this section, we build a reference crystal as a BCC (100) QS. 

Figure 2 shows the geometry for the BCC (100) QS. We take a “unit cell” in which the atoms 

below the QS are allowed to move under isotropic transformations. The atoms above the QS, can 

support both volume and surface-surface separation transformation. Nearest neighbors to AQS 

- - - - l+y 1 a 

are denoted R and R y . A typical R y (out of four) is R y = a(l,l, )-— o(l, 1,1) = —(1,1,7). I n 

order to use equation (4) we will compare the gradient terms along the z direction that is, we take 

7 

h r and n e both equal to z. We need to evaluate cos(0) = R y -z- ^ = . For the neighbors R , we 

‘ + C 

V 2 4 


use the previous expression with the sign changed (because they point opposite to the z 
direction) and with y= 1. Then the left hand side of equation (3) becomes (there are four 
contributions of each kind): 
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Figure 2. — The AQS in the quasi-free surface (100). The dark volume can only 
experience dilations controlled by x. The light volume can experience both 
dilations and separations controlled by -y. 


Equivalent Gradient 




P -l e -aR + 


1 + Zl 

2 ' 4 


\aR- p)R } 


D - 1 —CXRy 

e y 


V3 


where R - — a n x , and R v - a n x, — + — . 

2 0 r 0 -y. 


r 

4 


The left hand side of equation (3) can be derived more simply in a similar way: 

Equivalent Density = 4R p e " R + 4 R y p e 0X7 


(19) 


( 20 ) 


Application of the Method 

We use the method described in the Generalization of UBER section to predict surface energy 
formations of Fe (110) and (111). It is necessary to select the reference surface which has the 
largest density gradient, that is the second term in equation (1). Fortunately, this information is 
available from reference [2], With the advent in new ab initio techniques such information can be 
available for other surfaces of interest. 

In order to solve equations (3) and (4) we need evaluate the density and gradient terms for the 
real crystals of interest and compare them with expressions given in equations (19) and (20). 

1) Fe (1 10) surface energy formation. 

Iron’s ground state structure is BCC. Figure 3 shows a typical atom. When creating a surface 
(110) passing through that atom, two near neighbor atoms are lost. 
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Free surface 



Figure 3. — An atom in the real crystal loses two near neighbors when creating a (110) 
BCC surface. 


The contribution to the density from neighboring atoms is 

Density^ 110) = (21) 


with R 0 - ~^ a o ■> ar| d ao the equilibrium lattice parameter. 

To calculate the gradient term for BCC (110), we refer to figure 3. We notice that by symmetry 
the individual gradient terms from neighboring atoms in the (110) plane, when added together, 
cancel. Then the gradient term reduces to the contributions from the two remaining neighbors in 
the first layer below the free surface. By symmetry (they make the same angle with respect to the 
normal to the (110) plane) they both contribute to the sum on the right-hand side of equation (3) 

the same amount. One of those neighbors is at j? = -—(1,1,1). The direction of the nonnal to 


(110) is n = 


(HO) 

VI 



tenn becomes 


Since the two atoms contribute the same, the gradient 


Gradient (l 10 ) = — 2 



- P)K 


( 22 ) 
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Next we use equations (19), (20), (21) and (22) and substitute in (3) and (4). From that 
we find the values of y and x that satisfy the equations. Define 

A(x, y) = Equivalent Density - Density (1 10) ] 2 + [EquivaIent Gradient - Gradient (1 10)] 2 

Figure 4 shows a gray level plot of A(x,y) for the parameters of iron. From there a pair of roots is 
identified at x = 0.959 and y= 1.119. Substituting these values in equation (18) provides 
AE A qs = 2.49 J/m , which compares well with 2.247 J/m~ from first principles [2, 3], 

2) Fe (1 1 1) surface energy formation. 

Figure 5 shows a typical atom. When creating a surface (111) passing through that atom, four 
near neighbor atoms are lost. 



°-< U 

0.6 0.8 1.0 1.2 1.4 

x 

Figure 4. — Function to find the optimal values of x and y. 
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Real atom 

Figure 5. — An atom in the (111) Fe surface loses four near neighbors as compared with bulk. 


R 1 

The gradient term involves three atoms in the first layer below the surface with pq • h = — , and 


R 


one atom in the second layer below the surface with 


R 


R 


■ h = -1 . Then 


Gradient(\ 1 1) = — (3-j + l. !)(«/?„ — p)R£ 'e aR " =— 2 (aR 0 - p)R£ 'e a,t " 


(23) 


and 


Density{\ 11) = 4 R% e 0X0 (24) 

Next we use equations (19), (20), (23) and (24) and substitute in (3) and (4). From that 
we find the values of y and x that satisfy the equations. Define as in 1) 

A(x, y) = yj[ Equivalent Density - Densityi 1 1 1)] 2 + [ Equivalent Gradient - Gradient (1 1 1)] 2 

Figure 6 shows a gray level plot of A(x,y) for the parameters of iron. From there a pair of roots is 
identified at x = 1.07006 and y= 1.19797. Substituting these values in equation (18) provides 
AE A qs = 2.59J/m", which also compares well with 2.62 J/nT from first principles [2, 3]. 
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1.8 



0.6 0.7 0.8 0.9 1.0 1.1 1.2 


x 

Figure 6. — Function to find the optimal values of x and yfor the formation of the 
(111) surface. 

Table 1 shows results for Fe low index energy formation calculated by various methods. In 
addition, the table shows the ratios between (110) and (111) surface energies. This number is 
relevant because it has been used in the past as a gauge of the predictive power of a theory. It is 
known [8] that the ratios of surface energies are lower than expected from simple geometrical 
arguments based on surface roughness. This is due to the smoothing of the electron density 
parallel to the surface. Our results for the (111) to (110) energies are within 10% of the ab initio 
results. 


Table I. — Comparison of surface energies calculated from various method. 


J/m 2 

OECT [5] 

AB INITIO 
[2,3] 

MEC (This work) 

Fe(110) 

1.78 

2.25 

2.49 

Fe( 111) 

2.48 

2.62 

2.59 

Fe(\ 10) 
Fe( 1 1 1) 

1.39 

1.16 

1.04 
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Conclusions 


This paper introduced MECT, whose main features are the introduction of a gradient tenn in the 
density and the possibility for anisotropic defonnations of the equivalent lattice. This approach is 
of interest in the calculation of surface energies in particular, and more generally for low 
symmetry geometries in which the gradient of the density is significant and therefore makes an 
important contribution to the total energy. We have applied the method to calculate low index 
surface energy fonnation of Fe and found an improved agreement with ab initio results. 
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